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Abstract 

An exactly solvable many-particle quantum system is proposed by adding 
some nonhermitian but PT invariant interactions to the Calogero model. 
We have shown that such extended B^ Calogero model leads to a real spec- 
trum which obey generalised exclusion statistics. It is also found that the cor- 
responding exchange statistics parameter exhibit 'reflection symmetry' pro- 
vided the strength of a PT invariant interaction exceeds a critical value. 
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I. INTRODUCTION 



As is well known, in quantum mechanics one usually chooses a hermitian Hamiltonian to 
ensure real energy eigenvalues of the corresponding Schrodinger equation. However, quan- 
tum mechanical systems characterised by nonhermitian Hamiltonians also play a significant 
role in many contexts like absorption of incident particles in nuclear physics, localisation- 
delocalisation transitions in superconductors and in the description of the defraction of atoms 
by standing light waves Recently, theoretical investigations on different nonhermitian 

Hamiltonians have received a major boost due to the remarkable observation that many 
such systems, whenever they are invariant under combined parity (P) and time reversal 
(T) symmetry, give real energy eigenvalues This seems to suggest that the condition 
of hermiticity on a Hamiltonian can be replaced by the weaker condition of PT symmetry 
to ensure that the corresponding eigenvalues would be real ones. However, till now this 
is merely a conjecture supported by several examples 0-0. Moreover, in almost all of 



these examples, the Hamiltonians of only one particle in one space dimension have been 
considered. 

The aim of the present work is to test the validity of the above mentioned conjecture 
for the cases of some exactly solvable many-particle quantum mechanical systems in one 
dimension. For a Hamiltonian containing number of particles, the PT transformation is 
evidently given by 

i — > —i, Xj —Xj, Pj — > Pj (1.1) 

where j G [1, 2, • ■ ■ , A], and Xj {pj = —i^) denote the coordinate (momentum) operator of 
the j-th particle. So we want to find out some exactly solvable nonhermitian Hamiltonians 
which remain invariant under the PT transformation (|1 . 1|) and investigate whether such 
systems would lead to real spectra. 

In this context one may note that the well known A^r-i Calogero model, which contains 
A particles on a line and is described by the hermitian Hamiltonian 
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ha=-It.^+1^'t.^3++o e , ,2 ^ (1-2) 



represents an exactly solvable system [|T^]. This type of exactly solvable models with long- 
range interaction have attracted a lot of attention due to their close connection with diverse 
subjects like fractional statistics, random matrix theory, level statistics for disordered sys- 
tems etc. [Il7|-p4|. Recently an integrable extension of A^^i Calogero model is proposed by 
adding some momentum dependent interaction to the Hamiltonian ( |1.2| ) and this extension 
of Calogero model is also solved exactly to obtain the energy eigenvalues as well as eigen- 



functions |^ . It is found that these energy eigenvalues are real and bounded below, in spite 
of the fact that the added momentum dependent interaction given by 

(,3) 

2 ^^^Xj -Xk \dxj dxkj 

where 5 is a real parameter, is not a hermitian operator. 

However, we now interestingly observe that the complete Hamiltonian {Ha + Hp) of 
the above mentioned extended Calogero model is indeed invariant under PT transformation 
( 11.11 ). Thus we find here an example of PT invariant many-particle system, which is not 
only exactly solvable but also leads to completely real spectrum. To obtain more examples 
of this type, it may be noted that the Calogero model associated with root system can 
also be solved exactly p6| , |27|| . Therefore, it is natural to ask whether this Bn Calogero 
model can be generalised in a PT invariant way so that the newly constructed model would 
remain exactly solvable and yield completely real spectrum. In Sec. II we try to answer this 
question by proposing an appropriate PT invariant extension of the Bn Calogero model and 
solving such model exactly. In Sec. HI we study some salient features of this extended B^ 
Calogero model and show that its spectrum can be interpreted through generalised exclusion 



statistics |28| as proposed by Haldane. Sec. IV is the concluding section. 
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II. EXACT SPECTRUM OF A Bm TYPE CALOGERO MODEL 



Here we propose an extension of the well known Calogero model [2^,27] as 



1 N f)2 1 N n ^ ^ N 2 _i_ 2 



^ j=l ^-^j ^ j=l ^ j=l -^j j,k=i -^k) 

+ 5 y 1 A + 5 y 1 (x—-x —\ (2 1) 

where gi, g2, Si, 62 are some real coupling constants. In the special case Si = 62 = 0, the 
Hamiltonian ( |2.1D reproduces the original Calogero model: 

I N r.2 I N n ^ 1 Af ^2 I 2 

HB--\T.l, + l-'T-^ + 'iEi, + s.T.7^i^- (2.2) 

^ j=l ^-^j ^ j=l ^ j=l j,fc=i K-^j ^kJ 

0+k 



It may be observed that, though the Hamiltonian violates hermiticity prop- 

erty due to the presence of momentum dependent interactions like i^i SjLi ;^gf- + 
^2 Z^j^fc ^33^5" (^iaf" ~ -^^^af^)' remains invariant under the combined PT transformation 



( |1.1|) . Therefore, it should be interesting to solve this extended Bj^ Calogero model and see 
whether it gives real energy eigenvalues. 

In this context it may be noted that, both v4jv_i and Bj^ Calogero models have been 



solved recently by mapping them to a system of free harmonic oscillators For the 



purpose of solving the Hamiltonian (|2.1| ) through a similar procedure, we conjecture first 
that its ground state is given by a Laughlin type wave function 

= n (^i - ^kYixj + XkY n 4 e"^ (2.3) 

j,k = l 1 = 1 

p and a being two real non-negative parameters which are related to the coupling constants 

gi,g2,6i,52 as 

gi = a{a - 1) - 2a6i ; g2 = 2p{2p - I) - Ap52. (2.4) 



By solving the quadratic equations ( |2.4| ) of a and p, it is easy to see that these two parameters 
take real values provided the coupling constants in the Hamiltonian satisfy 



gi>-{\ + 5,f: g2>-{\ + 52f. (2.5) 



So, in this article we shall consider the PT invariant Hamiltonian (|2.1|) only for the range 
of parameters compatible with the conditions ( p.5|) . Now if we use the expression ( p.3|) 
for a similarity transformation to the Hamiltonian ( p.l|) , it reduces to a simple 'effective 
Hamiltonian' like 

I N Q2 N Q ^ N ^ / Q d \ ^ ^ I d 

/ j=l (JXj OXj ^ j,^^ Xj Xj. \ OXj OXk ) Xj OXj 

Here 

Eg^ = '^+ coNd- + ujN{N - l)2p, (2.7) 

and the real valued parameters a, 2p are defined as a = a — 6i and 2p = 2p — 62 ■ 
It should be noted that, at the limit 51 = ^2 = 0, Hb (p.6|) reduces to 

^B = ^ + 2p ^ 2 _ 2 i^ja ~ ^Z!— (2-8) 

'^j^lOXj OXj j i^^j^Xj X^ \ OXj OXk / j^i Xj OXj 



where the parameters p and a are related with the coupling constants gi,g2 of the B 



N 



Calogero model (|2.2|) as 

g, = aia-l); g2 = 2p{2p-l). (2.9) 



So, by using the expression ( |2.3|) for a similarity transformation on the Calogero model 
( p.2| ) and applying the relations ( p.9|) , one can get the corresponding effective Hamiltonian 
( p.8| ). It has been found earlier that this effective Hamiltonian ( |2.8|) has a complete set 



of polynomial eigenf unctions I^Sj. Such polynomial eigenfunctions, which are completely 
symmetric with respect to variables x'^, are called as generalised Laguerre polynomials. Since 
these symmetric polynomial eigenfunctions are evidently free from any singularity, there 
exists an one-to-one correspondence between the nonsingular eigenfunctions of the usual B^ 
Calogero Hamiltonian and the generalised Laguerre polynomials. Now it is very interesting 



to observe that, after a trivial substitution of coupling constants given by 2p — > 2p and 



a ^ a, the effective Hamiltonian ( |2.6| ) for our extended B]\i Calogero model coincides with 
the effective Hamiltonian ( |2.6| ) for Calogero model. Therefore, the generalised Laguerre 
polynomials again represent a complete set of polynomial eigenfunctions for the effective 
Hamiltonian ( |2.6| ). Consequently, there exists an one-to-one correspondence between these 
polynomials and the eigenfunctions of the extended Bn Calogero Hamiltonian (|2.1|) . 

/^From the above discussion it is clear that, by following Ref. ||35|, one can directly solve 



the eigenvalue problem for T-Cb (|2.6| ) through generalised Laguerre polynomials. At present, 
however, we want to follow a different approach which attempts to solve this eigenvalue 
problem by mapping Ti^ to a set of decoupled harmonic oscillators. To this end, we notice 
that this effective Hamiltonian ( |2.6| ) may be expressed a.s Hb = + uS^ , where the B^ 
type Lassalle operator and the Euler operator are given by pSl 



2 dxj .^^^^ Xj \ dxj dxk J Xj dxj 

S^ = ^.,^. (2,10) 

It is easy to see that the above defined operators satisfy the simple commutation relation: 
[S''^, 5"^] = —25*^. By taking advantage of such commutation relation, we perform further 
similarity transformations to T-Cb and reduce it to the Hamiltonian corresponding to free 
oscillators as 

H,.ee = S-'HbS = -- E + ^^UY./. - — (2.11) 



where S = Soe^^^^ and Sq = e^^^ e'^'^ 



N s2 



Due to similarity transformations in ( p.6| ) and ( |2.11| ), one may naively think that the 



eigenfunctions of the extended B^ Calogero model ( p . 1| ) can be obtained from those of 

the free oscillators as: ^/'ni.na, ,njv = ilJgrSo[Hn^{xi)Hn2{x2) ■ ■ ■ Hnj^{xn)) , where Hn^{xj) 

denotes the Hermite polynomial of order rij. However, similar to the case of Bn Calogero 
model ipOl,^, the action of «So leads to a singularity unless all rijS are chosen to be even 



integers and Hn^{xi)Hn^{x2) ■ ■ ■ Hn,^{xn) is symmetrised with respect to all coordinates. 
Therefore, nonsingular eigenfunctions of Tis ( |2.1|) can be obtained from the eigenfunctions 
of free oscillators as 

= l/jgrSo^^A^ {H2nAxi)H2n2{x2) ■ ■ ■ H2nr,{Xn)) ) , (2.12) 

where A+ completely symmetrises all coordinates and thus projects the distinguishable 
many-particle wave functions to the bosonic part of the Hilbert space. Eigenvalues of TYb 
corresponding to eigenfunctions ( p.l2| ) are given by 



En,,n2,...n^ = + 2cu ^ n,- = — + uNd + uN{N - l)2p + 2uY, (2.13) 
where the excitation numbers rijS are non-negative integers obeying bosonic selection rule 

Since a and p are real parameters, the energy eigenvalues ( p.l3|) are also real ones. 



Thus we interestingly find that the PT invariant nonhermitian Hamiltonian ( |2.1| ) generates 
real spectrum. However, within the above mentioned approach, it is not clear whether the 
corresponding eigenfunctions form a complete set. It is also evident that, apart form a 
constant shift for all energy levels, the spectrum ( |2.13| ) coincides with the spectrum of N 
number of free bosonic oscillators with frequency 2u. As the term 2ujY,jLin'j in eqn.( p.l3|) 



is always non-negative, -Eni,n2,.-.n]v '^^^ t)e less than Egr. Consequently, ipgr ( p.3|) indeed 



represents the ground state wave function of Ti.B (|2.1| ) with energy Egr. 

Finally we want to point out an important difference between the ground state energies 
of B]\f Calogero model ( p.2|) and its PT invariant extension ( |2.1|) . At the limit 5i = 62 = 0, 



eqn.( p.7|) reproduces the ground state energy of -Bat Calogero model as 

Egr = '^ +u;Na + u;N{N -l)2p, (2.14) 

and the form of corresponding ground state eigenfunction is given by ( |2.3| ). Since o and p 
must be positive to ensure the nonsingularity of ground state eigenfunction (|2.3|), the ground 
state energy (|2.14| ) of the Calogero model is always a positive quantity. However, as 



evident from eqn.( p.7|) , the ground state energy of the extended Bn Calogero model is 
dominated by the A^^ order term for large values of A^. The coefficient of this A^^ order 
term, i.e. 2p, is related to the coupling constants ^2, 92 through a quadratic equation given 
by 

^2 = (2p)'-2p-52(l + <^2). (2.15) 



By using ( p.l5|) and the second relation of (|2.4| ), it is easy to see that the parameter 2p 



becomes negative (even though 2p remains positive) if the coupling constants 62, g2 are 
chosen within the range: ^2 > 0, > (72 > —^2(1 + '^2)- Therefore, for large values of A^, the 
ground state energy of extended Bn Calogero model ( |2.1| ) will also be a negative quantity 
within the above mentioned range of two coupling constants. 

III. SOME PROPERTIES OF EXTENDED Bn CALOGERO MODEL 
A. Connection with fractional statistics 



Generalised exclusion statistics (GES) introduced by Haldane is believed to play 
an important role in the edge excitations of fractional quantum Hall effect. Such exclusion 
statistics can be realised microscopically in usual Calogero models with hermitian Hamiltoni- 



ans 



22|-p^. The GES parameter for these Calogero models is a measure of 'level repulsion' 



of the quantum numbers generalising the Pauli exclusion principle 1^^. So the partition 



function and various thermodynamical quantities for these Calogero models can be derived 
as a function of the corresponding GES parameters. Now for exploring GES in the case of 
our PT invariant Bn type Calogero model ( p.l|) , we observe that eqn.( p.l3l) can be rewritten 
exactly in the form of energy spectrum for A^ free oscillators as 

ON ^ . . 

where uj = 2uj and 

rij = nj + 2pj + h^2a -Sp-l) (3.2) 



are quasi-excitation numbers. However, from eqn.( p.2|) it is evident that such quasi- 
excitation numbers are no longer integers and they satisfy a modified selection rule: 
fij+i — fij > 2p . Thus the minimum difference between two consecutive fijS is given by 

2p = 2p-62, (3.3) 

(here we assume 2p > 62)- As a consequence the spectrum of extended Calogero model 
( ^31) satisfy GES with parameter 2p. It is interesting to notice that this exclusion statistics 
parameter depends only two coupling constants 62 and (72, which control the strength of 
long-range interactions in the Hamiltonian (|2.1| ). Moreover, eqn.( p.l5D describe a parabolic 



curve on the (^2, 92) plane for any fixed value of 2p. Consequently, all points on such a 
parabolic curve, representing extended Bj^ Calogero models associated with different values 
of 62 and g2, yield the same exclusion statistics parameter. 

It may be observed that the eigenfunctions ( p.3|) and (|2.12|) pick up a phase factor 



( — 1)^'' under the exchange of any two particles. So the exchange statistics parameter for 
the extended B]\r Calogero model (^TTj) is given by 2p. It is clear from eqn.(|3.3|) that, for the 
case 62 ^ 0, the exchange statistics parameter for the extended Bj^ Calogero model differs 
from the corresponding exclusion statistics parameter. Though the the exchange statistics 
parameter is not directly related to the spectrum or thermodynamics of the system, it 
has some interesting features. For example, in absence of confining harmonic potential, 
the scattering phase shift for multi-particle scattering generally depends on the exchange 
statistics parameter p4[. Moreover, quite similar to the case of usual Calogero models ||24|| , 



the exchange statistics parameter for the extended Bj^ Calogero model fixes the boundary 
conditions on the wave functions (|2.3|) and ( p.l2|) at the limits Xi Xj. 



B. Reflection symmetry of the exchange statistics parameter 

It is well known that, for fixed values of all coupling constants, the exchange statistics 
parameter of the B^^ Calogero model (|2.2|) may take two distinct values. These two distinct 
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values of the exchange statistics parameter are also related through a 'reflection symmetry' 
given by: 2p ^ 1 — 2p. For exploring such reflection symmetry in the case of extended i?Ar 
Calogero model we note that the second relation in ( |2.4| ) can be easily solved to obtain 
two solutions of the exchange statistics parameter as 

1 



2p = -(1 + 2^2) ± -V(l + 2^2)2 + ^92. (3.4) 

Due to the existence of these two solutions, it appears that reflection symmetry is also 
present in the case of extended Bn Calogero model (p.l| ). However, we have already noticed 
that only positive solutions of the parameter 2p lead to physically acceptable nonsingular 
eigenfunctions of Tis- So, the reflection symmetry can exist in the case of extended 
Calogero model provided both solutions of 2p in eqn. (|3.4|) take non-negative values. It is 
easy to see that both of these solution will be non-negative only within a parameter range 
given by (1 + 282) > and §2 < 0. Thus we curiously flnd that a kind of 'phase transition' 
occurs at the line 62 = — | on the (62,92) plane. For the case 62 < — |, the reflection 
symmetry of the exchange statistics parameter is lost for any possible value of 92- On 
the other hand, for the case 62 > —■^ the exchange statistics parameter shows a reflection 
symmetry: 2p — > 1 + 262 — 2/9 , if g2 is chosen within an interval — (| + 62)^ < 5'2 < 0. 



C. Relation with Calogero model 

We have seen in sec. II that, similar to the case of B^ Calogero model, the extended Bp^ 
Calogero model can also be mapped to a system of free harmonic oscillators through 
a similarity transformation. So it is natural to enquire whether the extended Bj^j Calogero 
model (|2.1|) is directly related to the Bn Calogero model (|2.2| ) through some similarity 
transformation. Investigating along this line, we flnd that 

1 1 /^' 1 Af ™2 _|_ ™2 

T-^HsT = - ^ E + E + I E i + ^2 E ^ > (3-5) 

3+k 

where 

10 



N , N 

.2 ^ ^2N% TT^^l 



r= n(^?-^D^n4\ (3.6) 

j.k = l 1=1 

and H'q denotes the Hamiltonian of Calogero model with 'renormahsed' couphng con- 
stants given by 

g[= gi + 5i{l + 5{), g'^ = 92 + H'^ + ^2) ■ (3.7) 
However it should be observed that, for any nonzero value of ^2, either F (|3.6|) or its 



inverse becomes singular at the limit Xj Xk- So the relation (|3.5|) can not be interpreted as 
a similarity transformation in the usual sense and it may lead to some strange consequences. 
For example, due to relation ( |3.5| ), one may expect that the Hamiltonians T-Cb and H'^ 
must share exactly same eigenvalues. However, we have already noticed in sec. II that the 
ground state energy of H'^ is always a positive quantity, while the ground state energy of 
Hb ( p.l|) can be a negative quantity provided the exclusion statistics parameter 2p takes 
a negative value. To explain this rather unexpected result, we first recall that 2p will 
be negative if the coupling constants of TYb ( |2.1|) are chosen within the range: ^2 > and 
> (/2 > — 62(1 + 62). So, from eqn.( |3.7|) one finds that the renormalised coupling constant g2 
must be a positive quantity in this case. Consequently, the corresponding exclusion statistics 
parameter 2p', which is obtained by solving the quadratic equation g2 = 2p'{2p' — 1), has 
one positive and one negative solution. The positive solution of 2p' leads to a ground state 
for H'q with positive eigenvalue. On the other hand, for large values of A^, the negative 
solution of 2p' yields a ground state for H'^ with negative eigenvalue. However, one usually 
throws away this negative solution of 2p', since the corresponding ground state eigenfunction 
ijjg^{xi,X2, ■ ■ ■ ,xn) becomes singular at the limit xj — > Xk- Nevertheless, as will be shown 
below, this unphysical eigenfunction of H'^ enables us to construct a physical eigenfunction 
of Hb- E)ue to the relation ( |3.5|) , one finds that 



'l/jgr{Xi,X2, - ■ ■ ,Xn) = T'l{)'g^{xi,X2, - ■ ■ ,Xn) , (3.8) 

where ipgri^Xi, X2, • ■ • , xat) represents the ground state eigenfunction of Ti.B- It can be easily 

checked that the zeros of F are sufficient to cancel all singularities of "0' ( ) at 
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the limit xj Xk ■ Therefore, the ground state eigenfunction ijjgr{xi,X2, • ■ • , xn) ( |3.8| ) is no 
longer singular at this limit. Thus we curiously find that a singular eigenfunction of H'^ gen- 
erates a nonsingular eigenfunction oi TCb through the transformation (|3.8|). Consequently, 



the ground state energy ofHs ( |2.1| ) coincides with the negative eigenvalue of H'^ associated 
with its unphysical eigenfunction '^^^.(xi, X2, • ■ • ,Xn). In a similar way one can show that 
all singular eigenfunctions of H'^, which represent the excited states over the unphysical 
ground state eigenfunction ipL-i 

), generate nonsingular excited states of Tin 

( p.l|) through a transformation by F. 



IV. CONCLUSION 

Here we observe that a recently considered nonhermitian variant of ^at-i Calogero Hamil- 
tonian with real spectrum remains invariant under the PT transformation. Being encouraged 
by this observation, we propose a new many-particle quantum system ( |2.1| ) by adding some 
nonhermitian but PT invariant interactions to the Bn Calogero model. Such PT invariant 
interactions depend on both coordinate and momentum variables of the particles. By us- 
ing appropriate similarity transformations, we are able to map this extended Sat Calogero 
model to a set of free harmonic oscillators and solve this model exactly. It turns out that 
this many-particle system with nonhermitian Hamiltonian yields a real spectrum. This fact 
supports the conjecture that the condition of hermiticity on a Hamiltonian can be replaced 
by the weaker condition of PT symmetry to ensure that the corresponding eigenvalues would 
be real ones. It is also found that the spectrum of extended Calogero model obeys a 
selection rule which leads to generalised exclusion statistics (GES). However, the extended 
Bn Calogero model also possess some remarkable properties which are absent in the case of 
usual Bn Calogero model. For example, we curiously find that the GES parameter for this 
extended Bn Calogero model differs from the corresponding exchange statistics parameter. 
Moreover a 'reflection symmetry' of the exchange statistics parameter, which is known to 
exist for B^ Calogero model, can be found in the case of extended model only if the strength 
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of a PT invariant interaction is chosen above a critical value. As a future study, it will be 
interesting to search for other exactly solvable many-particle systems with nonhermitian 
but PT invariant Hamiltonians and investigate the nature of their eigenvalues as well as 
eigenfunctions. 
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